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Abstract 

The small-x behavior of the singlet contributions to the polarized structure function 
92{x, Q"^) is calculated in the double-logarithmic approximation of perturbative QCD. The 
dominant contribution is due to the gluons which, in contrast to the unpolarized case, 
mix with the fermions also in the small-x domain. We find a power-like growth in 1/x in 
the odd-signature parts of the amplitude with the same power as in the singlet function 
gi{x, Q^) at X <C 1. 
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1 Introduction 



The investigation of the structure functions gi and g2 provides the basis for the theoretical 
description of polarization effects in deep inelastic lepton nucleon scattering. The Q^-evolution 
of the spin dependent function gi is well known: at x ~ 0(1) one can use the original Altarelli- 
Parisi equations , whereas in the region of very small x the double logarithmic approximation 
p|, 1^ has to be used and predicts a stronger growth in 1/x than the Altarelli-Parisi equations. 
The situation with the other polarized structure function, g2, is more complicated. Compared 
to gi, several new features appear Q]. 

First, there exists a twist-3 contribution to g2, g^\ which is not suppressed even for large 
1^. Part of the function g2 may be obtained from the known structure function gi{x, Q"^) with 
the help of the Wandzura-Wilczek relation p, and from another integral relation 0] which 
relates the twist-3 piece of g2 to the (m^/Q^-suppressed) twist-3 corrections to gi. However 
these relations are not enough to restore the whole spin dependent function g2{x, Q^). This can 
be seen from the explicit one loop calculation of g2{x,Q^) on a quark target P, ^ Q and also 
from the small-x behaviour of the nonsinglet component of g2{x,Q'^), calculated in the double 
logarithmic approximation (DLA) in [Q. 

Another distinct feature and difference from the standard leading-twist evolution equations 
is the fact that the number of operators which contribute to the twist-3 component of g2ix, Q"^), 
g^\ is not fixed but increases with the moment index n (see, e.g. 0)" Nevertheless, in 
the two limits: 

i) n — *• cxo (where n denote the (x)-moments) 

ii) number of colours —>■ oo 

it has been shown that the quark-gluon operators decouple from the quark operator evolution 
equation |]13[, and the asymptotic behaviour of g2{x,Q^) in the region 1 — x ^ 1 has been 
derived 0. 

Next, in the region of small x g2 has to be decomposed into pieces with different signatures. 
Whereas unpolarized structure functions and the polarized structure function gi have a definite 
signature (e.g., gi is the energy discontinuity of an odd-signature scattering amplitude), g2 
contains both even and odd signature [|r^, |T^, and an important part of the small-x analysis 
of g2 is the separation of the two signature structures. Because of the different contents of 
signature, gi and g2 are expected to have a different behaviour near a; = (i.e. different powers 
of see, for example, |18|). 

Finally, in the region of small x another new feature arises. It is well known that the 
strong ordering of transverse momenta is violated 0, |TB|, |T^, Instead of kf^ S> one 
has kf^ 2> -^^kf_^^^ (with Xi -C Xi^i). As a result, double log contributions of the form 

In^ l/xY appear which cannot be summed up in the framework of the conventional logQ 
evolution. This feature holds for both gi and g2-, and for gi it has been discussed in detail 

^Thc one loop function 32(2;, Q^) on the gluon target was calculated in 

^See also |l^for a recent study of the twist-3 (52'^^) evolution in the N^—f 00 limit and for a possibility 
to study the scale dependence of 32(2;, Q^) beyond the Nc ^ 00 limit. 
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in 0, ^. In particular, the infrared evolution equation (IREE) [|TB| has been used to sum up 
all the leading pQCD double logarithms (DL) of the form a" In l/xln'^Q^ with m < n and 
k + m = 2n — 1. InQ] similar techniques have been used to investigate the small-x behaviour 
of the polarized structure function so far only the nonsinglet case has been analyzed. In 
this paper we present the study of the singlet part of g2, and we find its asymptotic behaviour 
at x — >■ 0. In analogy with our study of the nonsinglet part of g2 we consider the scattering 
amplitude of a virtual photon an a quark. We will not (yet) address the question how our 
calculations can be used for the polarized structure function g2 of the proton, i.e. how to set 
up the initial conditions coming from the confinement region. 

In 0] we have shown that, for the even signature part of the nonsinglet g2, the only double 
logarithmic contributions come from the ladder-type Feynman diagrams with t-channel inter- 
mediate states consisting of two quarks only. For the odd-signature part of the amplitude one 
has to consider also non-ladder graphs where an additional 'soft' t-channel gluon "embraces" a 
part of the ladder diagram, in which the transverse momenta of all partons (gluons and quarks) 
are larger than the transverse momentum of this "embracing" soft gluon. The DL contribution 



of such a graph can be summed up by the IREE, using the method proposed in |T8|. In we 
have analyzed all one and two loop diagrams, contributing to the nonsinglet component of g2, 
and we extracted the double logarithms. In the singlet part of g2, whose analysis is the goal 
of the present paper, there are only two new 2-loop diagrams, shown in Fig.l. Fig. la has the 
same ladder structure but with the two-gluon t-channel intermediate state while the crossed 
box diagram of Fig. lb, as usual, looses the logarithm. The analysis of these graphs, together 
with the results of |^ will enable us to find the small-x behavior of the singlet g2 structure 
function. 

A comment is in place about the contribution of t-channel states with more than two gluons. 
For unpolarized structure functions at very small x, diagrams with any number of t-channel 
gluons contribute to the small-x behavior of the scattering amplitude ~ l/x (modulo powers of 
In 1/x). For the polarized case (both gi and (72) the antisymmetric tensor structure (see below) 
reqires, for the two gluon state, the polarization of the t-channel gluons to be different from 
each other; as a result, the amlitude behaves as ~ x° (modulo powers of In 1/x). Consequently, 
even at small x, quark and gluon t-chanel states mix, quite in contrast to the unpolarized 
case. However, also in the polarized case, there exist t-channel states of gluons which lead to a 
small-x behavior ~ 1/x. The simplest one PT| consists of three gluons, has positive C-parity, 



and is dual to the Odderon solution discussed in |^ . This three gluon contribution, however, 
is beyond the DL approximation used in this paper. First, the coupling of this state to the 
virtual photon is suppressed by a few extra powers of which are not accompanied by the 
maximal number of logarithms. Next, the contribution of this state to g2 is suppressed by 
an extra power of l/Q^, needed to compensate the dimension of the transverse momenta ku, 
k2t (of the gluons) which saturate the antisymmetric tensor £t^'"^l^^ and reflecting the orbital 
momentum in the three gluon state. For these reasons, contributions from diagrams with more 
than two gluons will not be considered in this paper. 

The outline of the paper is as follows. In Sect. 2 we consider the kinematics of the process and 
outline our strategy. In Sect. 3 we calculate the DL small-x contribution of the simplest singlet 
diagrams shown in Fig.l. It turns out that the DL asymptotics of Fig.l is proportional to the 
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singlet structure function gi^XjQ"^). Based on this general property we consider in Sections 3 
and 4 the small-x behaviour of the spin dependent function g2{x, Q^). Sections 5 and 6 contain 
the evolution equations and their solutions. In a concluding section we present a brief summary. 



2 Kinematics and definitions 



Let us begin with a brief summary of previous results and with an outline of our future strat- 



egy. The spin dependent part of the hadronic tensor iW^^ of the DIS lepton-quark scattering 
amplitude has the form 



where M,p^ and denote the mass, the four momentum and the spin (pseudo)vector of 
the target, respectively. As usual, (sp) = and = —1, and the spin-dependent structure 
functions are defined as the discontinuities of Ti and T2: 

9i,2 = - ^ Im Ti,2 . (2) 

Throughout this paper we work in the Feynman gauge and use the Sudakov representation of 
the momenta of quarks and gluons: 

ki = -aiQ + Pip + kti . (3) 

The photon four momentum is written as g = g' — xp, and p and q' are the two light-like 
reference vectors. The Jacobian can be written as: 

S 

d^ki = — daid(3i(fkti . (4) 

Here S = 2q'p denotes the center of mass energy, and we assume that the quark target mass is 
small: <^ Q"^ = -g^. So one may neglect the value of this mass M everywhere, except for the 
mass term in the target density matrix. This mass term is needed in order to obtain a nonzero 
spin trace proportional to the first power of M (see Eq.(P) . Finally, we recall from [^] that, 
in order to obtain in (|I]) a nonzero contribution to g2, we need a transverse polarization of the 
target (s^ = s^t): longitudinal polarization of the target = p'^/m contributes to the structure 
function gi only, while the transverse (s^ = Sf^) polarization gives the sum g± = gi + g2- 

Let us briefly recapitulate our previous analysis of the polarized structure function. Consider 
a generic ladder diagram (Fig. 2) and begin at the bottom. Since the spin dependent part of 
the hadronic tensor Wfj^iy is antisymmetric (proportional to the e-tensor), in the target density 
matrix we have to start from the 75S part: 

^ (p + m)(l -75s) . (5) 



^ Hence we may put ~ p 



/2 
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Figure 1: The lowest order Feynman diagrams which contain a gluon t-channel state 



When calculating the trace, this 75-matrix gives the antisymmetric tensor e^^*^"/^. As we have 
said before, for the present analysis we chose transverse polarization of the target quark. 

Moving upwards inside the ladder diagrams, we encounter t-channel intermediate states 
which consist of an quark-antiquark pair or of two gluons. For the quark t-channel state, we 
have found in [Q two different polarization structures which we will denote by 's' and 'k': 

's' , (6) 

'k' {2sk) ■ -f^k , (7) 

where fc^ is the t-channel quark momentum. At the target we begin with the 's' structure, but 
at any gluon rung above we have both possibilities: either to continue with 's' or start a 'k 
branch'. The latter continues with 'k' structures only, until at the upper end of the diagram it 
ends with the coupling to the photon. There is no return from a 'k branch' to an 's branch'. 
Morover, the 'k branch' contains only odd signature, and its evolution is the same as in gi, 
whereas the 's branch' contains both even and odd signature and has slightly modified evolution 
equations. Nonladder structures cancel for even signature branches, whereas for odd signature 
they lead to an additional term in the evolution equation. 

What is new in the present singlet case, is the two-gluon intermediate state. We first 
note that this intermediate state must have the same antisymmetric structure (with respect 
to the gluon polarizations /i' and z/') as the e polarization tensor. Thus we cannot assign 
to both gluons the same longitudinal polarizations which would lead to the small-x behavior 
~ 1/x. The best one can do for retaining the largest power of 1/x is to assign the longitudinal 
polarization {e^> oc q'^,) to one of the two gluons and the transverse polarization {e^i = e^'t) 
to the other gluon. Therefore, as mentioned before, in contrast to the unpolarized structure 
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function Q"^) ~ the prediction for the small-x behaviour of the polarized distributions 
goes as 

fl'i,2(a;, Q^) ~ const • In In as In^ , (8) 

where the function F denotes the leading (DL) corrections to the asymptotics of the Born 
singlet spin dependent structure function: (71 2(0;, Q^) ~ const as x ^ 0. Speaking in terms of 
the complex angular momentum j plane, both the two t-channel quark and two gluon exchange 
lead to the rightmost singularity at j = u; = 0, and even at very small x for the singlet spin 
dependent structure functions we have to account for the transitions between the quark and 
gluon t channel states. Continuing with the spin structure of the t-channel two gluon state we 
have to recall that the longitudinal gluon polarization e^i is contained in the term q'^ip'^n / {q'p') 
in the spin part of gluon propagator 



ip'q') 



where the vector g^, {p'^n) is contracted with the spin structure at the lower (upper) end of the 
i-channel gluon. The transverse polarization is contained in f?^/^//. At first sight there are two 
antisymmetric spin structures with one longitudinal gluon: 

(a) e^'^^Psip'p (10) 

and 

(b) s^p'^ - p'^s^ . 

However, the last structure has unnatural parity {s^ is a pseudovector) and never contributes 
to our amplitude Eq.(|I|). So it is enough to consider the first structure, Eq.(|TUp, only. Since 
this structure is the same as in the case of (71, we can make use the results of ||^. 

It will be the central task of this paper to compute the evolution kernels of the different 
t-channel states and the transitions from one to another. Altogether, we have three different 
states: the two spin structures of quark-antiquark states, (6) and (7), and the two gluon state. 
As usual, the DL-contribution comes from the kinematical region of 

1 > /9i > /92 > ••• > a; 

and 

. . . < ai < a2 < . . . < 1 . (11) 

The contour of the ai loop integration will be closed around the pole of the propagator of the 
i-th s-channel (horizontal line) particle. Thus all the s-channel particles lie on mass-shell, and 
the t-channel propagators can be approximated by 1/kf ~ /fcf^, as ctj/^jS -C k1^. So the i-th loop 
integral takes the form 

C — ^^-j^ ■ numerator . (12) 



2tx Pi k 



The logarithmic (3i integration is written explicitly, while, in order to obtain the second log- 
arithm (from the integral over rffc^J, we need to select the term proportional to /c| in the 
numerator, coming from the spin part of the diagram. As discussed before, the gluon— *^gluon 
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Figure 2: A generic Feynman diagram for the scattering of a virtual photon on a quark target 



case we take from 0; the evolution of the 's branch' and 'k branch' quark structures were 
obtained in the transition from 's branch' to 'k branch' is also contained in [Q, and the 
transition from 'k branch' to 's branch' vanishes. Collecting these results in terms of a matrix 
whose columns (rows) are labelled by the spins structure of the t-channel states below (above) 
the transition, we have for the numerator in (|T2l): 



99 kg sg 
gk kk sk 
gs ks ss 



( ^ 


-1 ? 


2 


1 




-1 1 



(13) 



Here the labels ((?, fc, s) refer to the spin structures of the gluon, the 'k' structure of the quark, 
and the 's' structure of the quark, resp. The 2x2 block matrix in the upper left corner coincides 
with gl and has been taken from |^ ; the 2x2 block matrix in the lower left corner, on the other 
hand, is the same as in the nonsinglet g2 case and has been derived in 0|. What needs still 
to be computed are the matrix elements denoted by '?': in the two following sections we will 
show that both elements vanish. 

The color factors C are the same as in the unpolarized scattering: for the gluon— i>gluon it 
is equal to Ca = Nc {Nc = 3 is the number of colours), for the quark —* gluon and quark 
quark splitting we have Cp = {N^ — l)/2Nc, and for the gluon —>■ quark transition the color 
factor is 2Tp = rip {np is the number of light quark flavours). Combining these color factors 
with the elements of the numerators matrix in eq.(13), we have all ingredients to the infrared 
evolution equations. 
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3 Singlet two-loop diagram 



We start with the amphtude Fig. la and use the transverse vector st^ The spin dependent lower 
part of Fig. la takes the form: 



Bu'u' = — Tr 



7,,(p + m)75S7^.(p -k + m)]= -218^^'^'''^ s^kp , (14) 



where in order to obtain the leading (in 1/x) contribution we have to keep the longitudinal 
component of the gluon momentum ~ jSp'^. 

Next, to simplify the calculation of the trace of the upper quark loop we make use of gauge 
invariance. The amplitude which corresponds to the upper block and includes both the 

box (Fig. la) and crossed box (Fig. lb) diagrams, satisfies the properties 

-q, = (15) 

and 

U^^."k," = U^,Ah + Pp' -aq'),. = 0. (16) 

Based on the gluon gauge invariance ([1^) one may replace the momentum p'^„ in the longitu- 
dinal gluon polarization vector (i.e. spin part of the gluon propagator (H)) by —kt^"/j3; the 
contribution of the term —aq'^„ does not contain the double logarithm and may be neglected 
within the kinematical domain ([TTD . We contract the tensor in (14) with the gluon numerators 
(9) and make use of this replacement. The polarization of the t-channel gluons is then described 
by the tensor structure: 

e^"kt^" - ktf,ne^u {ks^) _l {ks^) , , 

where e^j, is the two dimensional antisymmetric tensor, and the polarization vector has the 
form e^'i = ej^'ia^a- this way the expression (|T3p can be rewritten as 

- 2^(3e^"''"-^p'^si ^^^^^ ^'"'"'''p'J, ■ (17) 

Since the transverse spin vector s^^ corresponds to the spin-flip amplitude from the photon 
side we choose one transverse and one longitudinal photon polarization vector, i?^ = Ej^ and 
= EI = — i=(g' + xp'). With the help of the photon gauge invariance ( pTS]) the last vector 



can be replaced by Ej, = 2xp' l\fQ^. Thus, finally, in the case of Fig. la we have to calculate 
the trace 

Tr'^ = '^^^e^''^''-^p^J,{ill&{q^l^^^^^ 
Using the identity 

Iv^lii = —^^livXa^xlhla + (the tcrm symmetric under /i ^ z/) (19) 
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we obtain 



4:{ks^ 
{p'q') 



{p'{i - A;))Tr l-f^q'iE^^iq + i)E^ 



- {q'{i-k))TT -ir^p'iE\\q + i)E^ 



(20) 



Putting in E^ = 2xp'/^, it is easy to see that the largest (and the only DL) contribution 
(where q' is multiplied by the momentum k, and p' by momentum i'; note that ^ and 
^ Pi) comes from the second term and reads: 



4:{ks^ 
{p'q') 



{q'k){2p'i)TT k^llg^^ 



(21) 



As the s-channel quark with momentum (k — i) is on mass shell, the product {q'k){2p'i) / {q'p') 
—j3ka£S ~ (fc — i)j. So after the azimuthal angular integration we get 



that is 



4(A;s^)(A;-£)2£^ 



Tr'^ ~ 2kf£fTi j^E^E^s' 



IE' 



'8k 



(22) 



where the first factor is exactly the structure we are looking for in the hadronic tensor (|T]), and 
the last factor Sk^if is needed to keep the leading logarithms in the dkf and dif integrations, 
as it was discussed at the end of Sect. 2. 



Finally, let us consider the trace corresponding to the crossed box diagram of Fig. lb: 



Tr" 



2i 



{ks^ 

{p'q'] 



-e'^^^f'p'J^TT 



lu{q + i - k)E\i - k)j^iE^{q + . 



(23) 



The product of the three 7- matrices {iE^{q + 1)) may be decomposed as 

iE^{q + i) = cs^5 + dsi5ls. (24) 

The last term has unnatural parity and may be neglected, while the vector cg = ^TT[iE^{q + 
i)'ys]- Using this fact and the identity ([T9|) we get 

(ks^) 



Tr" 



{p'q') 



Tr 



p'iE^{q + i) -Tr -f5q'{q + i-k)E\i-k) 



Tr 



q'iE^{q + i) -Tr -f^p'^q + i - k)E\\i - k) 



(25) 



The second trace in the first term on the rhs of (|2^) vanishes as it contains only one transverse 
vector (i — k)t, while the last term in ( ^3]) is zero since in the second trace Ej^ oc p'^. Therefore 
we neglect the crossed box contribution (after the trick based on the gauge invariance was 
applied) and the whole DL result, given by the ladder-type Fig.l diagrams, reads: 



-Cg ■ 2cf 



Ols\' 

2tt) 



dPk 

/3k 



(26) 
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(here is the electric charge of the quark). In (26), the hmits of integration follow from (11) 
and the discussion after (11). 

We have to note that, within the DL approximation, this result for g± = Qi + g2 (0) 
coincides with the lowest order singlet function g^. In other words, to this order of accuracy 
we obtain the singlet (^) ~ 0- This is the consequence of the fact that in the lower part of 
Fig.l graph the Born structure function (7^ = 0. For the one loop approximation (order as) it 
is known P that at small x the one loop result is given by: 

92ix,Q ) ~ ^^xln—- 



i.e. is much smaller than 



2tt m?x 



gi[x,Q ) ~ In. 



2n rn?x 

So it is only at the two loop level (or above) that the nonsinglet structure function g'^''^' becomes 
comparable to the value of g^'^' at x -C 1. Thus one expects a non-zero DL contribution to the 
small- a: behavior of the singlet g2 component at earliest at two loops. Our result (26), however, 
shows that at the two loop level it still vanishes. We therefore have to turn to three or more 
loops. On the other hand, we know already from the argument given in section 2 that in the 
DL approximation the gluon contributions both to g± = gi + g2 and to gi are driven by the 
same spin structure (|T0|). Therefore we may expect that the distributions (/I and gl will have 
the same asymptotic small-x behavior. In the following section we will show that this is indeed 
the case. 



4 Higher order contributions 

We now turn to higher order corrections to Figs, la and lb. As a generic diagram, we consider 
Fig. 2. Beginning with the quark target at the bottom, we first have the same evolution as in the 
nonsinglet case. After the first s-channel gluon the initial quark density matrix |(]3+m)(l— 75s) 
gives 

^Il7a(p + m)(-75s)7a = m75S, (27) 

where we omit the terms which do not depend on the quark spin vector s^. In this way we 
obtain the structure 75S. At the next step this spin structure, together with the t-channel quark 
propagators k , produces the second structure {2sk) ■ 75A;: 

r75sr = (A;"2)75S - (2sA;") -75^". (28) 

Consequently, starting from two s-channel gluons, the quark density matrix in the upper {k') 
cell of the block contains the two structures 75s and {2s-^k')'-^^k' . So our diagram splits into 
two branches, and we have to follow each branch seperately. The latter one, 75A;', is equivalent 
to the density matrix of longitudinally polarized quark with the momentum k'^. So this spin 
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dependent part of the heavy photon-quark k' amphtude is described by the photon-quark 
structure function giipf*^) 

T,Al*q) = jj^ Qa-k'p- Tiix', Q') (29) 

with Im Tl = —27Tgi{'y*q) and x' = Q'^/2{k'q). The integration over the azimuthal direction 
of the vector k[ (/ d(p{2s^k^)k'^ = kfs-^) replaces the tensor e^^'^^gc^fc^ by e'^''"^qasj^ and gives 
the factor k'^ needed to save the DL structure of the dk'^ integral As a result, the branch in 
which originates from the 'k' quark structure, 2{s^k')'yr,k' , is proportional to g^. So we can 
make use of our results for gf: we know the transitions from the gluon t-channel to the quark 
t-channel (and vice versa), and we also know the latter never contains the 's' spin structure. 
In other words, we have shown that that in (13) the upper right element, sg, vanishes. 

Following the other branch, which starts from the 's' structure, 755^, at the upper end of 
the lower block, we note that it does not lead to any double logarithm at all. To show this 
we again make use of the gauge invariance, that is the condition analogous to (p!6D but with 
respect to the lower block: B^iyik^i = 0. Based on this condition we replace the vector q' by 
—kf^'/ak- After this the polarization of the only t-channel gluon state, which may give the 
DLogs is described by the pure transverse tensor T^',^/ oc e^/^,/ ■ {ks-^)/a. On the other hand, 
the 755-*- quark structure produces the tensor B^/^' = —Aie^''^'"^s-^{k' — k)f3 orthogonal to ej^/^'- 
As a result, in DL we have no transition from the 's' structure to gluon states: in the lower left 
element, gs, vanishes. 

In summary, the only DL singlet contributions start from the 75A;' quark structure, and the 
gluon evolution is the same as in g^. It therefore may be written as the convolution of the 
nonsinglet structure function g2 with the singlet function g^. Note that we have to use not 
the whole DL nonsinglet function g2'^', but only its part g^ which corresponds to the structure 
75/c'. This part can be easily extracted from ref.jQ, while the DL singlet function gf needed in 
the convolution can be taken from 

This value of ujg is considerably larger than the value of corresponding "intercept" of the 



nonsinglet function g2''^', where Un.s. = t^o = y C^Cpcts)/'^ (the ratio of the two values, lOs/ujo is 
about 2.6). Therefore the asymptotic behaviour of the spin dependent singlet structure function 
g2{x,Q'^) is almost completely driven by the DL behavior of the singlet gl{x,Q'^) distribution. 

Let us return to the signature question. As we have mentioned in section 2, the analysis of 
the nonsinglet case in has shown that the 's' branch of the quark state contains contributions 
to both signatures, wheras the 'k' branch has negative signature only. Since in our DL the two 
gluon state communicates only with the 'k' branch, it connects to negative signature only. The 
only place where even signature can appear is the 's' structure of the quark state which is the 
same as in the nonsinglet function g2'^'- 

The negative signature in the gluon t-channel state can be seen rather directly. Indeed, as 
it was discussed in [Q, one obtains the negative signature term keeping the photon (or gluon) 
momentum in the trace over the quark loop. In the quark gluon transition the trace reads 



XrS? = Tr 



75k'ju'{k'-k)'j^> , (30) 
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where the gluon momentum k gives the only nonzero contribution. The same can be seen in 
the gluon to quark transition. Recall our two loop calculation of Sect. 3. Here the leading DL 
contributions comes from the second term of the rhs of Eq. (pOD . In this term we have to keep 
the product (q'k), which includes the gluon momentum k and the photon momentum q in the 
trace (see (|2T|) ) in order to obtain the double logarithms. 

Let us finally comment on a point which, at first sight, might look dangerous. Namely, 
starting from the spin structure 

i^-^'k^P, ^ (31) 



(pk) (pk) 

in the lower nonsinglet part of the Fig. 2, where the gluon momentum ka plays the role of the 
photon momentum q in Eq.(|l]), it looks as if one might get a much larger small-x behavior, 
g2 ~ l/x, than obtained in the previous discussions. In comparison with (|I^ this term 2(5-*"^) • 
ie^^'^^kaPp — k1 ■ ie^^'^^ Sa ■ Pp is enhanced by the factor 1/ (3k as the longitudinal component of 
gluon momentum ka — Pkp'a- However an extra factor l/jSk destroys the logarithmic structure 
of the d/Sk/Pk integral. Even more, we will show that such a contribution vanishes. Indeed, 
the integration in Fig. 2 over the rapidity of the upper horizontal quark line of the lower block 
dPk' / Pk' may be written in terms of the x' = k"^ /2{pk) ~ k^/2{pk) variable. This last quark 
with momentum {k' — k) is on mass shell. So, keeping the transverse momenta fc^ and kf fixed, 
we have 

It-' _ i-|2 1 1.' _ t.|2 

I \t I \t 



k' 



akS 2p'k 
and therefore dPk' / Pk' = dx'/x'. 

On the other hand, we invoke the argument which leads to the Burkhardt-Cottingham sum 
rule [p^ 



dx'g2{x',Q'^) = . (32) 

This sum rule can be derived from an unsubtracted disperison relation for the photon - proton 
forward scattering amplitude. For the very special hilicity structure which contributes to g2, this 
amplitude leads to a cross section which falls with energy faster than 1/s; therefore, the integral 
/ ds'g2{s') is finite, and from the large-s behavior of the dispersion relation one finds (32). The 
same large-l/x condition applies to our starting function g2'^- (shown in |Q), therefore we 
conclude that 

^ g^^x',^) T^Y = / dx'g{x',e) = (33) 



x' ' (pk) 

Thus the term which looked so dangerous, in fact, does not contribute to double logarithmic 
structure function g2. 



5 Infrared evolution equations 

With results of the previous sections we now turn to the infrared evolution equations (IREE). 
We are interested in the scattering of a virtual photon on a transversely polarized quark target: 
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= Ti + T2, and we will concentrate on the odd signature part. In order to formulate 
our coupled evolution equations, we first have to generalize to different target structures. As 
outlined in section 2, inside our diagrams we have to consider two different spin structures 
in the quark t-channel state (||),(0), labelled as and g^, and the gluon spin structure ([T0|). 
Correspondingly we introduce the three component vector: 

(Thg) \ 

T = n^qk) . (34) 

V tUqs) J 

For our scattering of a virtual photon on transversely polarized quark target we will be interested 
in the g^-component. We write T as a Mellin transform: 



T = / TT- K a^)Ri^,y), (35) 




where R{uj^y) is a three -component vector, defined on analogy with (0), and uj = j denotes 
angular momentum. The signature factor ^{uj) = has the form: 

(36) 



and 

y = \n{^\. (37) 




The structure function g']_ is obtained from the discontinuity (H), and we have to take into 
account both DL-contributions and zvr-terms. Before we can write down the TREE, we have 
to introduce further auxiliary quantities: scattering amplitudes Fij (with i,j = g,qk,qs) which 
describe the scattering of quark or gluon with spin structure z on a target with structure j 
(loosely speaking, these scattering amplitudes are obtained if, in Fig. 2, we remove the coupling 
to the photon at the upper end). We write them as Mellin transforms, i.e. the Fij are partial 
waves and depend upon angular momentum u. By Fq we denote the 3x3 matrix composed of 
the Fij (the subscript refers to the color singlet quantum number): 



^ ^99 ^kg Fgg ^ 

Fgk -^fcfc Fg}^ 
\ Fgs Fjfg Fgg j 



(38) 



(it is the exact analogue to Fq in 0] and to /d ■* in p|). 

We now turn to the evolution equations. The TREE (illustrated in Fig. 3) describes the 
variation of the amplitude with respect to the infrared cutoff /x^: 

= i^ + l)R- (39) 



^^J? 

This differential operator stands on the the left-hand side of the TREE for T, which is illus- 
trated in Fig. 3. The right-hand side of the IREE is derived from the observation that the 
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Figure 3: Structure of the infrared evolution equations (IREE) 



dependence upon the cutoff /i resides in the intermediate state with lowest virtuahty (Fig. 3): 
the /^-derivatives of the amphtudes are equal to R times quark or gluon scattering amplitudes 
Fij with the external legs having transverse momenta close to /i. In terms of the matrix Fq, 
the evolution equation for the vector R becomes: 

The matrix Fq satisfies a nonlinear evolution equation (Fig.3): 

Foiu) = ^ Mo - GoFsH + FoM^. (41) 

Here the matrix Mq contains the evolution kernels (13) which we have computed in the previous 
section or collected from earlier studies: 

/ ACa -2Tp \ 
Mo = \ 2Cf Cf (42) 
V -Cf Cf J 

As we have discussed before, for the gluon state (upper left corner) the infrared evolution 
is described by the same splitting function as the singlet structure function gf. For the two 
quark system in the lower right two by two block matrix we can use the results of 0]. The 
remainig matrix elements describe transitions from quark to gluon states, they are taken from 
the previous section. The second term on the rhs of (^) belongs to the gluon bremsstrahlung 
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diagrams. The matrix Go has the form: 





Go = \ Cf I . (43) 



\ Cf 



In analogy to the matrix Fq which carries color zero in the t-channel we define the matrix Fg of 
color octet amplitudes. Due to the antisymmetric color structure the elements of Fg are even 
signature amplitudes. They satisfy evolution equations |T8| similar to (p]): 



Fs = ^ Ms + 1^ Fsiou) + -\- Fsiouf. (44) 
uj Stc^uj auj Stt^coi 

The color factor Ca in front of the second term on the rhs is the analogue of the matrix Gq in 
(PI). The difference between Ca in (|4^) and Gq in (|4l|) is due to the fact that for the positive 
signature amplitude Fg the sum of the two bremsstrahlungs diagrams is independent of the 
type of the incoming partons, and the matrix of color factors Gg becomes Ca times the unit 
matrix. Finally, the matrix Mg is the analogue of the matrix Mq in (RTI), but for color octet 
quantum numbers in the t-channel instead of color singlet. It reads: 

2Ca —Tf 

Mg = I Ca -l/2iV, I . (45) 
l/2N^ 




6 Solution of the evolution equation 

The solution for g]_, that is our vector T (eq.(34,35)), is obtained by solving Eq. (^) for Fg, 
then Eq.(PD for Fq and Eq.(?D) for R (the latter makes use of the Born approximation, Rb-, 
as initial condition for i?), and finally inserting R into (pS)). The final answer for the three 
component vector T will have the form: 

J 2m \xj \P / a; — Fo/Stt^ 




Rb = 2el (47) 



where 



Let us go through these steps in somewhat more detail. 

We begin with the equation for Fg. We first diagonalize the Born term, i.e. the matrix 
Mg. Thanks to the vanishing of the off-diagonal elements (Mg)^^ = (Mg)^^ = the two 
(largest) eigenvalues coincide with the analogue eigenvalues for the DL evolution of the 
singlet function gl. They are given by: 

= ± \^{2Ca + l/2Nf - ACaTf (48) 



15 



and 



A 



is) 



2Nr, 



(49) 



Let e^^\e^ ^ and e*^*-* denote the eigenvectors of Ms- 

I 1 



,{+) 



and 



X 



{+) 



A(+) - M- 



11 



,(+) 



12 



{-) 



22 



,(~) 



21 



AW - A(+) 



AW 
A(«) - A(-) 



(50) 



(51) 



The last eigenvector corresponds to the pure quark eigenvalue A = Ag and reads: 






1 



(52) 



Denote by = {e^^\e^ \e^'^'^) the matrix composed of these eigenvectors. For the diagonal- 
ization we also need its inverse E7^: 



1 



1 



1 - 





1 

5(+) J 



(53) 



where 



b^+) = _ and b^~^ = y^+^ x^-^ - y^'\ (54) 

Then we can diagonalize the matrix Mg: 

Mg = EsMsE^' , (55) 

where Mg = diag(A8 , Ag , Ag''"'). Consequently, eq.(P^ becomes diagonal if we transform to Fg: 

Fg = EsFsEs' . (56) 

Using the ansatz Eg = Ncq"^ lnM*(ci;) for the nonlinear Riccati equation (43), one finds second 
order linear differential equations for the -u*. Their solutions are given by parabolic cylinder 
functions. As a results we find for the components of Fg: 



^8 = N^g' ^\^e^"''D,Xz) 
where Dp denotes the parabolic cylinder function with 



(57) 



(58) 
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and 

UJ 



z = — , 000= jN.gySTT^ . (59) 

With this solution for Fg we return to the evolution equation(^TD for Fq which is solved by 
the (matrix-valued) expression: 



l-7T7^^o+ , ,, GoFs . (60) 



We will need the right-most singularity in the uj plane of the matrix Fq. Similar to the case of gi, 
this singularity is due to the vanishing of the square root in (|60|) , i.e. we need to determine the 
zeroes of the eigenvalues of the matrix under the square root. The diagonalization of this matrix 
is done in exactly the same way as for Mg: again, the matrixelements ^sg^ and 's/c' are zero, 
and the two largest eigenvalues are the same as for the singlet function gi. The diagonalization 
is done through the matrix Eq which consists of the corresponding three eigenvectors e^^^ , Cq ^ 
and eo"*^ (analogous to (50) -(52)). As it was discussed in 0, the largest eigenvalue A'^+^ is not 
too far from the value obtained in the pure gluonic case, neglecting the quark contribution. 
The accurate values have to be found from a numerical computation of the parabolic cylinder 
functions. For the case oi np = 4 we found in the rightmost singularity ojs = 3.45a;o (the 
pure gluonic case would have given 3.66 instead of 3.45, or even 4 if one neglects the non-ladder 
DL contribution). With ujq from (|59|) and = 0.18 we find for the leading term in ( ^6]) 



Us = 1.01 . (61) 

Two other singularities correspond (for = 4) to z^''^ = 1.81 and z^ = 1.39. 
Having found the matrix Fq, we return to (40) and find the matrix R: 




i? = (^^J R{u) (62) 

where the matrix valued function R{uj) has to be determined from the initial conditions of the 
evolution equation, from the energy dependence of R at the point y = 0, i.e. at Q"^ = yU^. At 
this point the matrix R satisfies another evolution equation (which has been discussed in detail 
in ^ (after eq.(3.14)), and its solution has the form: 

Riu) = Eo 1-—E,'Rb. (63) 

u — Fq/Stt'' 

Retaining in Fq only the leading upper component we find for the behavior of R near the square 
root branch point at u = Us'. 

~ T^W^l \-(^T'^ (64) 
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with = 0.29, x(^) = 0.43, = -0.052, and y^-^ = -1.40 Without the non-ladder 
contribution (i.e. putting (Fg = 0)) we would have x'-^-' = 0.28 and y^'^'^ = —0.039. 

Finally, to obtain the structure function g2 we have to subtract the known function gl 
{92 = 9± — 9i), and to add the positive signature part {—g^/2) given by the same pure ladder 
DL evolution as in the nonsinglet case 

9^ = 91 + 92= f ^ (-Y ^uj) (%]'^ E, l-—E,'Rs-9'/2 (65) 



7 Conclusions 



In this paper we have studied the small-a; behavior of the singlet polarized structure function 
g2, using the double logarithmic approximation. Our analysis is based upon ladder-type Feyn- 
man diagrams, but in order to make our DLA study complete, we had to take into account 
also the bremsstrahlung of soft gluons. We have shown that, at small x, the leading contri- 
bution to (72(^5 Q^) is given by the convolution of the nonsinglet spin dependent distribution 
(corresponding to the quark density matrix {s'^k)'y:^k) with the singlet function g^. As a novel 
feature, g2 contains pieces with different signature, both for the nonsinglet and for the singlet 
case. At small x, we therefore have to decompose the Feyman diagrams, in order to seperate 
even and odd signature. The small-x asymptotic behaviour of the structure function comes 
from the odd-signature part of the amplitude and is driven by the asymptotics of the singlet 
spin dependent function gf. At 1/x 00 we have found 

where the anomalous dimension 7 has the value 

7 = (67) 

and 

Us ~ 3.45^a,A^c/27r (68) 

(the latter result holds for = 4). In contrast to the unpolarized case, at small-x the 
gluons mix with the quarks. As in the singlet polarized structure function gf, the dominant 
contribution comes from the t-channel two gluon states. This is mainly due to the color charge 
of the gluon which is much larger than that of the quark. Numerically the value of ujg is 
rather large: Ug — 1.01 for = 0.18. This value of Ug is considerably larger than the value 
of corresponding "intercept" of the nonsinglet function 5'2'*'- ^^^'^ larger than the even 

signature part —g^/2 of the singlet g2, where ujn.s. = ^0 = \J {2CFCis)/T^ (the ratio of the two 
values, oJsIujq is about 2.6). Therefore the asymptotic behaviour of the spin dependent singlet 
structure function g^ix^Q"^) is almost completely driven by the DL behavior of the singlet 
gl{x,Q^) distribution. 
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From the formal point of view, the double logarithmic approximation (DL) that we have 
been using here (i.e.the neglect of nonlogarithmic contributions) can be justified only for the case 
of a very small QCD coupling cts ^ 1, in which case also cug is small. The experience from the 
LO and NLO BFKL calculations supports the belief that, whenever the leading approximation 
turns out to be large, one has to expect also large NLO corrections. Thus one should be 
careful in using the DL results for the singlet structure function for the numerical applications. 
Nevertheless, our findings indicate that the singlet spin dependent structure function g^l grows 
steeply at x — > 0. 
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